clear all; close all; clc;
R = 100;
C = 33e-6;
RC = R*C;
T = RC / 1000;
Tmod = 10*RC;
t = 0:T:Tmod;
%Проверка на иестовом воздействии
E = 0*(t<Tmod/4) + 3*(t>=Tmod/4);
U = nan(1, length(t));
dU = nan(1, length(t));

U(1) = E(1);
dU(1) = 0;
U2 = U;

for k = 2:length(t)
    %Вариант 1 - ток на резисторе определяет производную
    U(k) = U(k-1) + dU(k-1)*T;
    dU(k) = (E(k)-U(k))/RC;
    %Вариант 2 - результат замены производной на разность
    U2(k) = (RC*U2(k-1) + E(k)*T)/(RC + T);
end

figure;
plot(t, E, t, U, t, U2)
xlabel('t, sec');
ylabel('U, V');


